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Abstract 


Random  variables  T, ,T„,...,T  are,  in  this  paper, 

±  z  n 

associated  if  each  pair  of  non-decreasing  functions  F(T^,T2» •  •  •  »Tn)  » 
G(T^,T2» . . . ,Tn)  have  a  non-negative  covariance.  Association  holds 
in  cases  rangxng  from  Ti«T2*'**,Tn  independent  to  Ti»T2’*‘*,^n 
jointly  restricted  to  a  non-decreasing  curve.  Association  is 
preserved  under  the  standard  multivariate  operations  of  extracting 
subsets  and  po^xlng  independent  sets;  and  under  the  special  operation 
of  forming  sets  of  non -decree sing  functions.  Suitable  choices  of 
F,G  lead  to  various  inequalities  for  associated  random  variables. 

The  properties  of  association  are  studied  in  the  simple,  but 
representative,  case  that  Ti»T2*****^n  are  finitely  discrete. 

The  notion  of  association  is  useful  in  extending  the  domain 
of  validity  of  the  minimal  cut  lower  bound  for  the  reliability  of  a 
coherent  system  [4,1],  notably  (here)  to  the  case  of  repairable 
components  with  exponential  times  to  failure  and  exponential  times 

f 


to  repair. 


1 .  Introduction 


It  is  customary  to  think  that  two  random  variables  S,T  are 
associated  (in  the  sense  of  similar  behavior)  if  Cov[S,T]  ■»  EST-ES*ET 
is  non-negative.  The  same  two  random  variables  are  more  strongly 
associated  if  Cov  [  f (S) ,g(T) ]  >_  0  for  all  pairs  of  non-decreasing 
functions  f,g.  The  same  two  random  variables  are  still  more  strongly 
associated  if  Cov[F(S,T) ,G(S,T)]  >_  0  for  all  pairs  of  functions  F,G 
which  are  non-decreasing  in  each  argument. 

The  strongest  of  the  three  criteria  has  a  natural  multivariate 
generalization.  We  consider  the  definition 

(1.1)  Random  variables  Ti  *^2  ’  *  *  *  ,Tn  are  associated  if 

Cov[F(p,G(T)]  >  0,  where  X  -  <ti’t2 . Tn>  •  for  all 

pairs  F,G  of  non-decreasing  functions. 

Here,  and  later,  we  say  that  a  function  F  is  non-decreasing  if  it  is 
non-decreasing  in  each  argument,  i.e.  if  F(js)  <_  F(t)  whenever 
sA  <  t ^ i  i-l , ... ,n. 

This  paper  is  Intended  to  indicate  an  approach  to  the  organization 
of  the  properties  of  (1.1),  and  to  motivate  it  by  a  sample  application. 

We  study  (Section  2)  the  properties  of  (1.1)  in  the  case  that  the 
joint  probability  distribution  of  Ti»^2,’*',^n  is  finitely  discrete, 
i.e.  the  vector  X  takes  a  finite  number  of  values.  The  finitely  discrete 
case  illustrates  the  general  situation;  the  properties  of  (1.1)  obtained 
are  true  in  general ;  but  by  confining  our  attention  we  can  defer  discussion 
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of  some  equivalent  versions  of  (1.1)  which  are  of  interest  in  the 
general  theory. 

Our  application  (Section  3)  is  a  demonstration  that;  the  minimal 
cut  lower  bound  for  the  reliability  of  a  coherent  system  is  valid  in 
a  greater  variety  of  circumstances  than  previously  known.  This 
application  reflects  the  origin  of  our  Interest  in  (1.1),  and  is 
suitable  because  it  Involves  a  number  of  the  manipulations  that  our 
study  of  (1.1)  is  designed  to  facilitate.  We  do  not  attempt  a 
complete  discussion  of  this  subject.  Further  we  postpone  consideration 
of  certain  non-reliability  results,  e.g.  Kimball  (1951)  and  Robbins 
(1954),  which  can  be  related  to  (1.1). 

There  is  an  extensive  [see  Goodman-Kruskal  (1959)  for  a  recent 
bibliography]  and  continuing  literature  on  quantitative  measures  of 
association,  primarily  for  bivariate  distributions.  By  contrast,  our 
Interest  is  in  the  qualitative  identification,  and  corresponding 
properties,  of  classes  of  associated  distributions.* 

The  multivariate  dependence  described  by  (1.1)  is  one  of  several 
similar  notions.  When  referring  to  it  we  should  qualify  the  term 
"association",  but  for  the  present  it  is  simpler  to  let  it  mean  (1.1). 


*At  the  time  this  paper  was  being  typed,  our  attention  was  drawn  to  the 
manuscript  "Some  concepts  of  dependence"  by  E.  L.  Lehmann,  which  treats 
several  concepts  of  association  and  their  applications. 
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2.  Properties  of  association  for  finitely  discrete  and  binary  random 
variables 

For  finitely  discrete  Ti*T2** *  * »Tn»  functions  F(T)  have  finite 
expectations  EF(£).  The  covariances  considered  in  (1.1)  are  finite. 

2.1  General  properties 

Association  has  two  properties  desirable  in  any  classification  of 
multivariate  distributions: 

(Pj^)  Any  subset  of  a  set  of  associated,  random  variables  is  also  a 
set  of  associated  random  variables . 

(P2)  If  two  sets  of  associated  random  variables  are  independent  of 
one  another ,  then  their  union  is  a  set  of  associated  random 
variables. 


Property  P^  follows  Immediately  from  the  definition  of  association. 
P2  requires  a  short  proof,  which  can  be  made  by  considering  two 
independent.  Internally  associated  random  vectors  £  “  (S^ «s2 • • • • »sn) , 
T,  -  (T^Tj . Tm)  and  writing 


Cov[F(S,T) ,G(S,T) ] 


E£.TFG  "  E£,ZF 


E£,TG  " 


YxFG  -  W 


E_E_G 

&  X 


■  W*  -  h{V  ■  V  +  hlV  •  ■  W  •  W 

-  E^CoVjlF.G]  +  CoVg[ EjF , E^G ] , 


where  Ec  denotes  expectation  over  the  distribution  of  £.  E_  expectation 
over  the  distribution  of  T,  and  Ec  _  expectation  over  the  joint 
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distribution  of  £  and  £.  Eg  ^  -  EgE^  from  tha  independence  of  £ 
and  J£.  Since  Cov^tF^,J£)  »G(,g>,<2)  1  >_  0  for  each  fixed  £, 

EgCov,[F,G]  >_  0  follows.  Since  EjF(j^,T)  ,  EjGCA.Jf)  are  non-decreasing 
functions  in  £,  CoVg[ETF,E^G]  >_  0. 

Another  standard  multivariate  property,  valid  for  association,  is 

(Pj)  The  eet  consisting  of  a  single  random  variable  is  associated. 

To  show  P^  consider  a  random  variable  T  taking  values  t^  <  t2  <***<  t^. 
For  non-decreasing  F  and  G,  F(t^)  <_  F(t2)  <•••<  F(tk)  and 
G(t^)  G(t2)  <_•  •  •<_  G(tfc) .  It  follows  from  an  Inequality  for  similarly 

ordered  sequences  due  to  Chebyshev  [Hardy,  Llttlewood,  and  Pdlya  (1934), 

e 

12.17]  that  Cov[F(T) ,G(T) ]  >_  0. 

F3,  together  with  P? ,  implies  that  independent  random  variables  are 
associated.  The  example  of  independent  random  variables  represents  an 
"extreme"  of  association.  Another  extreme  is  represented  by  the  example 
of  random  variables  £  ■  (T^,T2 » ■  •  •  »T  )  taking  values  <•••< 

where  M,  <  &,  means  s^f^t^,  i"l,.«.»n,  and  s^  <  tA  for  some  1.  Then  fo 

F,G  non-decreasing  we  have  again  F(t^)£  F(t^^)  <_•  •  •<_  F(t and 
G(t(1))  <.  G^2^  <.*••<,  G(£*k)),  so  that  by  the  Chebyshev  inequality 
Cov[F(£),G(£)]  ^  0. 

Properties  P^,  P2  and  P^  permit  some  standard  multivariate 
manipulations  with  associated  random  variables.  The  applications  of 
association,  e.g.  to  reliability  theory,  are  often  founded  on  another 
manipulative  property. 
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(P^)  If  Ti»T2»***»Tn  are  associated,  then  any  set  of  non-decreasing 
functions  S^T)  ,S2(jT)  , . . .  ,Sm<jT)  are  associated. 

To  show  P^  we  observe  that  S^(Y)  ,S2(T)  , . . .  ,SmQ)  are  finitely  discrete; 
that  if  F,G  are  non-decreasing,  then  FQp,  G(S)  are  non-decreasing; 
and  that 

Cov^tF^S)  ,0(^)1  -  Cov^tFtS/Tn.G^T)}]  >_  0. 

2.2  A  working  definition  of  association 

In  the  finitely  discrete  case  there  is  one  alternative  definition 
of  association,  equivalent  to  (1.1),  that  is  particularly  useful. 

Theorem  2.1.  T2.,T2 . Tn  are  associated  if,  and  only  if, 

(2.2.1)  Cov[r<T),A(J)]  >0 

for  all  pairs  r,A  of  binary,  non-decreasing  functions. 

Proof.  Condition  (1.1)  clearly  implies  (2.2.1).  To  show  that  (2.2.1) 
implies  (1.1),  suppose  that  F(£)  takes  the  values  f^  <  f2  <•••< 
and  write  F(T)  -  ^  +  (f2-fi>r2<£)  +  •••+  (fk-fk_i>rkCE> »  where  I'i(X)  “  1 
if  F^T)  f^,  T^^T)  ■  0  otherwise.  Similarly  write 

G(T)  -  gx  +  (g2-g1)A2^T)  +•••+  (gJl-gJl_1)AJl(p .  Then 

Cov[F(T),G(T)]  -  Si_22J.2(fi"fi-l)<*j"*J-l>CovIritt>»Ajtt)l  -L  °* 

One  application  of  (2.2.1)  is  the  observation  that  for  two  binary 


random  variables  X,Y 


(2.2.2) 


Cov[X,Y]  >_  0  implies  X,Y  are  associated. 


This  result  follows  readily  from  (2.2.1)  since  the  only  possible  binary, 
non-decreasing  functions  r(X,Y)  are 

(r  -  X) 

(r  =  o)  <.  (r  -  xy)  <_  <  (r  ■  x  +  y  -  xy)  <  (r  e  i). 

(r  -  y> 

The  covariance  between  any  pair  of  binary  functions  r,A  such  that 
T  <_  A  is  automatically  non-negative,  and  the  remaining  pair  T  ■  X, 

A  -  Y  has  non-negative  covariance  by  hypothesis. 

The  utility  of  (2.2.1)  partially  rests  on  the  well-known  fact  that 
for  binary  random  variables  X,Y,  Cov[X,Y]  >_  0  is  equivalent  to 

(2.2.3)  P[X«1,Y«1]  •  P[X-0,Y-0]  >,  P[X-l,Y-0]  •  P[X-0,Y«1] . 


In  (2.2.1)  each  pair  of  binary,  non-decreasing  functions  r,A  defines 

a  partition  of  the  sample  space  of  the  vector  JT  into  the  four  disjoint 

regions  {r«j,A"k},  j ,k  -  0,1.  The  condition  (2.2.3),  applied  to  the 

pairs  of  binary  random  variables  r(,T)  ,A(T)  ,  gives  a  good  description 

of  the  joint  probability  distribution  of  associated  T, ,Tot...,T  in 

l  ^  n 

terms  of  the  geometry  of  the  sample  space. 

In  Section  2.1  we  mentioned  two  "extremes"  of  association.  Now 
(2.2.2)  and  (2.2.3)  permit  a  reasonably  complete  discussion  of  a  simple, 
non-extreme  example.  Let  a  vector  of  random  variables  T  take  values 


< 


*  £&  ■*(  i#k#  T-'-T 


****!*'. 


f 

1? 
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with  the  possible  additional  orderings  <.  or  i 

excluded.  In  this  case  we  can  define  X^(t^k)  “  j,  X2^t^k)  ■  k, 

j,k  ■  0,1,  and  have  X^(T),  X2Q)  non-decreasing.  We  can  also  write 
Ik 

T^  »  t^  when  X^  ■  j ,  X2  “  k  and  have  T^X^.Xj)  non-decreasing. 

It  follows  from  that  ' ’"^n  are  associated  if,  and  only 

if,  XX(T) ,  X2<T^  are  associated.  Then  from  (2.2.2)  and  (2.2.3) 

Xj^T),  X2(;i)  are  associated  if,  and  only  if,  P^Pqq  —  P01P10*  where 
pjk  -  P[X1(T)-j,X2(T)-k]  -  P[£-£,Jk]. 

2.3  Binary  random  variables 

Certain  properties  of  associated  random  variables  have  a  special 
intuitive  content  if  the  random  variables  are  binary.  The  most  useful 
of  these  is 

(BP^)  If  X1>X2,...,Xn  are  associated  binary  random  variablee3  then 
1-X^ , I-X2 , . . . , l_xn  are  associated  binary  random  variables. 

To  show  BP^,  observe  that  if  T  is  a  binary,  non-decreasing  function, 

then  the  dual  function  rD(jc)  ■  1  -  r  Ci-^)  ,  where  J.  -  jjc  ■  (1-x^ ,  l-Xj . 1-x^)  , 

is  also  binary  and  non-decreasing.  Let  Y  ■  ^  -  £.  Then 

covy[rgf),AQr)]  -  cov^[rD(jy  ,adqc)]  >_  o, 


and  the  property  follows. 

For  X1,X2,...,Xn  associated  and  r(£)  -  X^  A(X)  -X2X3...Xn,  (2.2.1) 

gives  E(X^X2* .  .X^)  >^E(X^)  •  E(X2X^ . . .Xr) .  Proceeding  inductively  we  find 

'■'rrt 

f 


n  u  p 

>_  n±-1  EX^*  Since  for  a  binary  random  variable  X, 

EX  -  P[X-1],  we  have 

(2.3.1)  PlXj-l.Xj-l . Xn-1)  >.  PtX^l]  •  PtX2-l]  •. . .  .P[Xn-l]  . 

From  (2.3.1)  and  BP^  follows 

(2.3.2)  PfX^O.X.-O . X  -0]  >  PtX.-O]  •  P[X,-0]  • . . .  *P[X  -0] . 

l  z  n  —  ±  z  n 

3.  A  sample  application  of  association,  the  minimal  cut  bound  in 
reliability  theory 

The  minimal  cut  lower  bound  for  the  reliability  of  a  coherent 
system  was  obtained  in  Esary-Proschan  (1963)  in  a  simple,  basic  case. 

In  this  section  we  show  that  the  bound  is  valid  in  some  of  the  more 
complex  situations  considered  in  reliability  theory.  The  derivation 
illustrates  the  application  of  a  formal  organization  of  the  properties 
of  association. 

3.1  Association  in  time  of  performance  processes 

In  the  setting  of  reliability  theory,  the  performance  process  of 
a  device  is  a  binary  stochastic  process  (X(t),t  ex),  t  C  [0,+<»), 

(e.g.  t  “  [0,+“)  or  t  “  {0,1,2,...}),  where  X(t)  »  1  if  the  device 
is  functioning  at  time  t  and  X(t)  ■  0  if  the  device  is  failed  at 
time  t. 

A  performance  process  {X(t),t  ex}  is  associated  in  time  if  for 
every  finite  set  of  times  {^^2*  •  •  ■  c  T»  the  binary  random 

. X(tk>  are  associated. 


variables  X(t^),X(t2) 
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The  basic  example  of  a  performance  process  that  is  associated 
in  time  arises  when  the  device  has  a  life  [ Esary-Marshall  (1964)],  i.e. 
when  P[X(s)  >_  X(t)]  ”  1  whenever  s  >  t.  In  essence  a  device  has  a  life 
if  once  failed,  it  stays  failed.  For  t^  <  tj  <•••<  t^  the  vector 
(X(t^) ,X(t2) , • • • ,X(t^))  can,  with  probability  one,  take  only  the  values 
(0,0,. ..,0)  <  (1,0,. ..,0)  <  (1,1,0, ... ,0)  <•••<  (1,1,. ..,1)  and  thus 

X(t^),X(t2) . X(t^)  are  associated  (cf.  the  cases  of  "extreme" 

association  in  Section  2). 

3-2  The  exponential-exponential  process 

Aside  from  the  cases  in  which  a  device  has  a  life,  perhaps  the  most 
frequently  studied  performance  process  is  the  alternating  renewal  process 
that  arises  when  a  device  falls,  then  revives  (e.g.  is  repaired  or 
replaced),  then  fails,  and  so  on,  where  the  time  S  from  revival  to 
failure  is  exponentially  distributed,  P[S  >  s]  ■  e  *8,  X  >  0,  s  >_  0,  and 
the  time  T  from  failure  to  revival  is  also  exponentially  distributed, 

P[T  >  t]  ■  e  Mt,  p  >0,  t  ^  0.  The  corresponding  (exponential-exponential) 
performance  process  (X(t),  t  e  [0,  +<*>)}  is  Markov  with  the  transition  law, 
for  s  <  t, 

P[X(t)-l|X(s)-l]  -  (A+p)_1{v+X  exp[-(A+y)  (t-s)  ]  } 

(3.2.1) 

P[X(t)-0|X(s)-0]  -  (A+u)-1{A  +  u  exp[-(A+p)  (t-s)  ]  }  • 

[See  e.g.  Cox  (1962)  or  Barlow-Pro schan  (1965)]. 

Theorem  3.1.  For  any  initial  distribution ,  i.e.  P[X(0)-1],  the 
exponential- exponential  performance  process  is  associated  in  time. 
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Proof.  We  consider  t^,t2» • • • >1^  such  that  0  <  t^  <  t2  <•••<  t^, 
and  let  ”  X(tj),  J  -  l,...,k.  We  show  that  Xj 9X2 9 • • • »X^  flre 
associated  by  induction.  alone  is  associated  by  P^.  Assume  that 

X^ 9 X2  9 • • •  jXj ar e  associat ed • 

Let  Pj  -  PfXj-llXj^-1]  and  qj  -  P[XJ-0|XJ_1-0]  .  From  (3.2.1), 

Pj  +  qj  i  (A+u)  ^(A+ij)  -  1.  Construct  binary  random  variables  Uj,Vj* 
Independent  of  X^»X2 . Xj-1*  with  the  joint  distribution 

P[Uj-1,Vj«1]  -  1  -  qJt  PtUj-l.Vj-O]  “  P-j  +  ~  1*  and  P[U  “O.V  -0]  -  1  -  p  . 

For  this  Joint  distribution  P[U,“1]  *  p  and  P[V  =0]  “  q  . 

J  J  J  J 


Figure  1 

Let  ej (Xj-i’^j »vj>  “  Xj-lUj  +  (1-Xj-i)Vj-  We  can  re8ard  xj  as 
generated  by  the  branching  process  shown  in  Figure  1  in  the  sense  that 
X^.Xj, . . • »Xj_^ »6j  have  the  same  joint  probability  distribution  as 

Xx , X2 9***9 X j  1  * X j  *  a ince  e . g . 
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P [ 9 j “1  »Xj _^“1  where  jjt  ■  (Xj_2»^j«3>'  •  • 

-  P[UJ-l,Xj_1-l%STLJ  -  PtUj-1]  •  PlXj.i-lJSTSl 

-  P[Xj“l  |Xj_1“l]  •  PCXj^-l^-xJ  -  P [Xj “1  ,Xj _1“1  ]  • 


Because 
associated, 
decreasing  in 


UjiVj’  Cov[Uj  ,Vj  ] 
Also  because  Uj  >_  Vj 
Xj_1 ,  U,,  and  Vj . 


> 


» 


0,  and  by  (2.2.2) 
the  binary  function 


Vvj  are 

0j  is  non- 


Now  X^,X2> . . .  ,Xj_^,Uj  ,Vj  are  associated  by  P2,  since 
X1*X2* * ‘ * ,Xj-l  are  associated,  Uj.V^  are  associated,  and  the  two 
groups  of  variables  are  Independent.  Then  X^.X^ . . . ,Xj_1,6j  are 
associated  by  P^,  since  these  variables  may  be  viewed  as  a  set  of  non¬ 
decreasing  functions  of  Xj^  ,X2 , . . .  ,X^_1  ,Uj  ,Vj  . 

Thus  X^,X2 . Xj_^,Xj  are  associated,  because  their  joint  distribution 

agrees  with  that  of  X^.X^. . .  ,Xj_^,0j  . 

The  same  method  can  be  used  to  show  association  in  time  for  the 
geometric-geometric  performance  process  (X(t) ,t*0,l, . . . },  i.e.  the 
alternating  renewal  process  where  time  to  failure  S  has  the  geometric 
distribution  P[S  >  s]  -  ps,  s  -  0,1,...,  and  time  to  revival  T  has  the 
geometric  distribution  P[T  >  t]  ■  qC ,  t  ■  0,1,...  .  Zn  this  case  the 
condition  p  +  q  >_  1  becomes  a  necessary  hypothesis,  but  the  result  is 
again  Independent  of  the  initial  distribution. 


■«'«*»■*  JW.V'  = 
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Other  performance  processes  of  Interest  In  reliability  theory 
are  also  associated  in  time,  e.g.  the  process  that  arises  from  scheduled 
replacement  of  a  device  whose  intrinsic  performance  process  is  associated 
in  time,  and  that  of  a  device  with  several  unloaded  or  lightly  loaded 
spares  (under  exponential  assumptions  about  times  to  failure  and  repair). 
For  our  present  purposes,  however,  we  want  to  focus  attention  on  the 
elementary  life  and  alternating  processes. 

3.3  Joint  performance  processes  generated  by  systems 

We  consider  systems  whose  performance  at  a  time  t  is  determined 
by  the  performance  of  their  components  at  the  same  time  t.  If 

QC(t),t  e  t},  where  £(t)  -  (X^(t) . Xn(t)),  is  the  joint  performance 

process  for  the  components  in  a  system,  the  performance  process  for  the 
system  is  {$£(t),t  e  x},  where  $jC(t)  =  $(£(t))  is  a  binary  function 
that  takes  the  value  one  If  the  system  is  functioning  at  time  t  and  takes 
the  value  zero  if  the  system  is  failed  at  time  t.  A  system  is  coherent 
[see  Barlow-Proschan  (1965)  for  a  survey  and  bibliography]  if  is 

non-decreasing  and  $(1,1, . . .  ,1)  ”  1,  $(0,0 . 0)  »  0.  The  function  $ 

is  called  the  structure  function  of  the  system. 

We  write  {$>X(t)  ,t  e  x},  where  $X(t)  -  ($xX(t)  ,$2&(t), . . . ,  $^(0  )  , 
for  the  joint  performance  process  of  a  set  of  m  systems,  with  structure 

ri 

functions  $i»$2’ * ' * »  ^m*  formed  from  a  set  of  n  components  with  the 
joint  performance  process  {£(t),t  e  x}.  We  admit  the  case  in  which  a 
component  participates  in  several  systems. 
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The  joint  performance  process  {X(t),t  e  t}  for  a  set  of  n 
devices  is  associated  in  time  if  for  each  finite  set  of  times 
{ tf, tj  i  • «  •  i  tfc}  C  t,  the  binary  random  variables  in  the  array 


(3.3.1) 


x^i)  x1(t2)  ...  x1(tk) 

x2<ti)  x2(t2)  *2^^ 


xn(ti)  xn(t2) 


x  (t.  ) 
n  k 


are  associated. 


If  n  devices  perform  independently,  and  the  performance  process 
of  each  device  is  associated  in  time,  then  the  rows  in  (3.3.1)  are 
independent  and  each  row  is  associated.  From  P2  their  Joint  performance 
process  is  associated  in  time. 


Theorem  3.2.  If  the  joint  performance  process  QC(  t )  ,t  ex)  of  a 
set  of  components  is  associated  in  time ,  then  the  joint  performance 
process  { j>X(t)  ,t  c  t}  of  a  set  of  coherent  systems  formed  from  those 
components  is  associated  in  time. 

Proof.  For  (tj.tj . tfe)  c  t,  each  $tX(tj),  i-l,...,m, 

j  ”  l,...,k,  is  a  non-decreasing  function  of  variables  in  the  array 
X^(tj),  i  -  1 . .  j  -  l,...,k.  The  conclusion  follows  from  P^. 

Using  Theorem  3.2  one  can  show  that  the  joint  performance  process 
corresponding  to  the  multivariate  exponential  (life)  distribution  defined 
in  Marshall-Olkin  (1966)  is  associated  in  time. 
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3.4  The  minimal  cut  lower  bound  for  the  reliability  of  a  coherent 
system 

The  reliability  at  time  t,  F(t)  ,  of  a  device  with  the  performance 
process  {X(t),t  £  t>  Is  defined  as  F(t)  -  P[X(s)»l,  all  s  e  -r(t)], 
where  r(t)  -  tO  [0,t].  We  make  the  convention  that  the  sanple  functions 
of  any  performance  process  are  continuous  from  the  right  on  t.  The 
reliability  of  a  device  then  becomes  the  survival  probability  or  the  time 
T  until  the  first  failure  of  the  device  occurs,  i.e.  P[T  >  t]  =  F(t), 
t  j>  0. 

Every  coherent  system  has  a  finite  number  of  minimal  cuts 
'  *  *  *cm’  A  minimal  cut  is  a  minimal  set  of  components  such  that  if 
all  are  simultaneously  failed,  the  system  is  then  failed.  The  performance 
process  {$jC(t),t  ex}  of  the  system  can  be  related  to  the  performance 
processes  {X^(t),t  e  t},  i  -  l,...,n,  of  its  components,  through  the 
minimal  cuts,  by 

(3.4.1)  4>jC(t)  -  rTj„1n.JX(t),  where  n^t)  -  1  -  n±eC  (1-X^t)), 

for  all  t  e  t  [see  e.g.  Barlow-Proschan  (1965)].  The  performance 
process  {hjjC(t),t  e  t}  is  that  of  the  (coherent)  system  in  which  the 
components  in  the  minimal  cut  perform  in  parallel.  The  representation 
(3.4.1)  corresponds  to  connecting  these  minimal  cut  parallel  structures 
in  series. 


The  minimal  cut  lower  bound  was  previously  obtained  in  what  amounts 
to  the  case  of  components  that  perform  independently  and  have  lives.  The 
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following  theorem  extends  the  bound  to  the  case  in  which  the  joint 
performance  process  of  the  components  is  associated  in  time. 

Theorem  3.3.  Let  F(t)  be  the  reliability  at  time  t  of  a 
coherent  system.  Let  Fj(t),  j  -  l,...,m  be  the  reliability  at  time 
t  of  its  minimal  cut  parallel  structures .  If  the  components  in  the 
system  have  a  joint  performance  process  that  is  associated  in  timet 
then 

(3.4.2)  F(t)  _>.  TT^PjtO.  for  all  t  ^  0. 

Proof.  By  Theorem  (3.2)  the  joint  performance  process  of  the 
minimal  cut  parallel  structures  is  associated  in  time.  Let 
Sk  -  {slts2, . . . .s^}  c  t (t) .  From  (3.4.1)  and  (2.3.1) 

(3.4.3)  F[tiCU£)-l,  A«l,...,k]  "  PlTTj®!  nj£<s£>"l»  ,k) 

*  •  •  •  »®i  i  TTj— ^  . « >  »k] * 

Now  let  +  S,  where  S  Is  a  set  dense  In  x(t).  By  monotone  convergence 

(3.4.3)  becomes 

P[4>(X(s)-l,  all  s  e  S]  LTTjml  P[n^X(s)-l,  all  s  e  S] . 

Then  (3.4.2)  follows  from  the  right  continuity  of  sample  functions  on  r(t). 

3 . 5  The  minimal  cut  lower  bound  as  an  approximation  to  system  reliability, 
a  numerical  example 

When  the  minimal  cut  lower  bound  ia  applicable,  it  is  simpler  to  find 
than  the  actual  system  reliability,  because  the  bound  requires  only  the 
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calculatlon  (or  simulation)  of  the  reliabilities  of  systems  of 
components  In  parallel.  Viewed  as  an  approximation  to  system  reliability 
the  bound  Is  conservative,  as  is  usually  desired  in  practice.  The  bound 
can  be  a  very  good  approximation  to  actual  system  reliability. 

A  numerical  example  is  given  in  Esary-Proschan  (1963)  in  the  case 
of  components  that  perform  independently  and  have  lives.  The  approximation 
is  good  when  component  reliabilities  are  reasonably  high,  the  most  likely 
practical  case. 

We  want  to  give  another  numerical  example,  in  the  case  that  components 
perform  independently  according  to  identical  exponential-exponential 
processes.  We  consider  the  "two  out  of  three"  system,  i.e.  the  system 
which  functions  whenever  two  or  more  of  its  three  components  are  functioning. 
The  system  has  three  minimal  cuts  -  {c^.Cg},  C2  -  {c^c^},  C3  »  (c^c^, 

where  c^,  c2,  denote  the  components.  In  this  case  both  the  reliability 

of  the  system  and  the  reliability  of  the  minimal  cut  parallel  structures  can 
be  calculated.  For  the  "2  out  of  3"  system  the  reliability  at  time  t  is 

(3.5.1)  F(t)  -  (k2~k1)  1(k2  exp [-k^a]  -  k^  exp[-k2a]}, 

where  2^  -  (b+5)  -  (b^lOW-l)**,  2k2  -  (b+5)  +  O^+lOb+l)*8,  and  a  =  At, 
b  “  p/A  [Halperin  (1964) J.  For  the  minimal  cut  parallel  structures 
(which  may  be  considered  "1  out  of  2"  systems)  the  reliabilities  at  time  t 
F^(t)  -  F2 (t)  -  F^(t)  are  given  by  (3.5.1)  with  2k1  «  (b+3)  -  (b2+Cb+l)'2, 
2k2  ”  (b+3)  +  (b^+bb+l)*5  {Epstein-Hosf ord  (1960),  Halperin  (1964)].  It 
is  assumed  in  deriving  (3.5.1)  that  all  the  components  are  functioning  at 


time  0. 
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From  (3.5.1)  F(t)  and  F^(t),  ^(t),  F^(t)  depend  only  on 
a  «  t  /  A  l.e.  the  ratio  of  mission  duration  to  the  mean  time  '  1 

to  component  failure,  and  b  »  X  ^/p  i.e.  the  ratio  of  the  mean 

time  X  1  to  component  failure  to  the  mean  time  p  ^  to  component 
revival.  In  Table  1  we  compare  F  with  *^’2  *F3  over  a  range  of 

values  for  a,b  that  are  of  likely  practical  interest,  i.e.  t  >  X  ^ 
and  X  ^  >  >  p  ^ . 

The  case  of  a  "2  out  of  3"  system,  when  the  components  have  lives 
distributed  according  to  the  trivariate  exponential  distribution  of 
Marshall-Okln  (1966)  ,  furnishes  an  example  in  which  the  minimal  cut  lower 
bound  is  a  poor  approximation  to  system  reliability — because  there  is  a 
positive  probability  that  two  or  more  components  will  fail  simultaneously. 
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s 

1 

10 

100 

1000 

10,000 

.44 

.089 

.0536 

.05016 

.049824 

2  a 

.36 

.074 

.0520 

.05001 

.049809 

— 

.565 

.336 

.3048 

.30156 

.301230 

■i 

.510 

.323 

.3034 

.30141 

.301216 

10  a 

.682 

.565 

.5505 

.548976 

.5488281 

.650 

.559 

.5498 

.548910 

.5488215 

50  a 

.8982 

.88799 

.88703 

.886931 

.8869215 

.8948 

.88758 

.88698 

.886927 

.8869211 

in2 

.9449 

.94205 

.941793 

.9417673 

.9417648 

10  a 

.9439 

.94194 

.941782 

.9417662 

<>417647 

103  a 

.994053 

.994021 

.9940183 

.9940180 

.9940180 

.994041 

.994020 

.9940181 

.9940180 

.9940180 

10A  a 

.9940053 

.9994002 

.9994002 

.9994002 

.9994002 

.9940041 

. >994002 

.9994002 

.9994002 

.9994002 

Table  1 

Comparison  of  system  reliability  and  the  minimal  cut  lower  bound 
for  the  "2  out  of  3"  system. 


The  components  perform  Independently  with  exponential 
times  to  failure  and  exponential  times  to  revival. 

a  -  duration  of  mission/mean  time  to  failure 

b  ”  mean  time  ->  failure/mean  time  to  revival 

The  bound  is  tabled  .  w  the  corresponding  system 
reliability. 
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4 •  Remarks  on  some  conditions  for  bivariate  dependence 

As  we  observed  in  the  introduction,  for  two  finitely  discrete  random 
variables  S,T,  each  of  the  conditions  in  the  list 

(4.1)  CovtS.T]  _>  0 

(^.2)  Cov[  f  (S)  ,g  (T)  ]  >_  0,  all  non-decreasing  f,g 

(4.3)  S,T  associated 

implies  its  predecessor.  We  can  add  to  the  list  the  condition 

(4.4)  T  is  stochastically  increasing  in  S, 

i.e.  P [ T  >  t | S“Sj ]  is,  for  all  t,  non-decreasing  in  j,  where 

s^  <  S£  <•••<  s^  are  the  possible  values  of  S.  It  is  not  difficult  to 

show  that  (4.4)  implies  (4.3),  by  an  adaptation  of  the  proof  of 

We  have  seen  in  (2.2.2)  that  for  two  binary  random  variables  X,Y, 
(4.1)  implies  (4.3).  The  construction  used  in  the  proof  of  Theorem  3.1 
can  be  readily  modified  to  show  that  for  binary  X,Y,  (4.1)  implies 

(4.4) .  Thus  in  the  bivariate,  binary  case  conditions  (4.1)  through  (4.4) 
are  equivalent. 

For  general  finitely  discrete  random  variables  S,T  none  of  the 
conditions  are  equivalent.  It  is  easy  to  find  S,T  satisfying  (4.1)  but 
not  (4.2).  If  S,T  each  take  the  values  a^  <  Sj  <  a^  and  have  the 
joint  probability  distribution. 
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8/64 

0 

15/64 

0 

18/64 

0 

15/64 

0 

8/64 

S“fl^  S*&2  3=83 


then  S,T  satisfy  (4.2)  but  not  (4.3).  If  S,T  take  the  same  values 
with  the  joint  distribution. 


T-*3 

1/8 

0 

1/4 

CM 

CO 

1 

H 

0 

1/4 

0 

T'ai 

1/4 

0 

1/8 

S=a2 

S=a3 

then  S,T  satisfy  (4.3)  but  not  (4.4). 
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